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ABSTRACT 

0^ | We consider the X-ray luminosity difference between neutron star and black hole soft X-ray 

transients (NS and BH SXTs) in quiescence. The current observational data suggest that BH 

H ■ SXTs are significantly fainter than NS SXTs. The luminosities of quiescent BH SXTs are 

^ | consistent with the predictions of binary evolution models for the mass transfer rate if (1) 

accretion occurs via an ADAF in these systems and (2) the accreting compact objects have 
event horizons. The luminosities of quiescent NS SXTs are not consistent with the predictions of 
ADAF models when combined with binary evolution models, unless most of the mass accreted in 
the ADAF is prevented from reaching the neutron star surface. We consider the possibility that 

^vq , mass accretion is reduced in quiescent NS SXTs because of an efficient propeller and develop a 

' model of the propeller effect that accounts for the observed luminosities. We argue that modest 

winds from ADAFs are consistent with the observations, while strong winds are probably not. 

m ■ 

t-H ' Subject headings: X-ray: stars - binaries: close - accretion, accretion disks - black hole physics 

- stars: neutron - stars: magnetic fields 

ON ■ 

6 ■ 

' Soft X-ray Transients (SXTs) are compact binary systems in which a low-mass secondary (either a 

main-sequence star or a subgiant) transfers mass via Roche-lobe overflow onto a black hole (BH) or neutron 
star (NS) primary (see reviews by Tanaka & Lewin 1995; van Paradijs & McClintock 1995; White, Nagase 
& Parmar 1995). SXTs have highly variable luminosities. They spend most of their lifetimes in a low 
luminosity quiescent state, but occasionally undergo dramatic outbursts during which both the optical and 
X-ray emission increase by several orders of magnitude (e.g. Chen, Shrader & Livio 1997; Kuulkers 1998). 
NS SXT outbursts typically occur every 1-10 years and last for several weeks, while BH SXT outbursts are 
typically separated by 10-50 years (or perhaps longer) and last for several months (see Chen et al. 1997). 

A variety of observations (see e.g. Tanaka & Shibazaki 1996) indicate that, near the peak of an 
outburst, SXTs accrete matter via a standard thin disk (Shakura & Sunyaev 1973), so that there is 
little doubt that the accretion is radiatively efficient during this phase. The situation is more complex in 
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quiescence. The spectra of quiescent BH SXTs do not resemble that of a thin disk, and the accretion rates 
inferred from the observed X-ray luminosities disagree by orders of magnitude with the predictions of the 
standard disk instability model for quiescent disks (e.g. Lasota 1996). 

Narayan, McClintock & Yi (1996) and Narayan, Barret & McClintock (1997) showed that the 
observations of quiescent BH SXTs can be explained by a two-component accretion flow model consisting 
of an inner hot advection-dominated accretion flow (ADAF; Narayan & Yi 1994, 1995a; Abramowicz et 
al. 1995; Ichimaru 1977; see Narayan, Mahadevan & Quataert 1998 and Kato, Mineshige & Fukue 1998 
for reviews of ADAFs) surrounded by an outer thin disk. In the recent version of this model described in 
Narayan et al. (1997a), only the inner ADAF contributes to the observed optical, UV and X-ray emission of 
the system. The outer thin disk acts mainly as a reservoir which accumulates mass until the next outburst 
is triggered. In quiescence, the emission of the disk is primarily in the infrared; this radiation is hardly seen 
as it is swamped by the emission of the secondary (see Narayan et al. 1996, 1997a and Lasota, Narayan & 
Yi 1996 for details). 

A key feature of the Narayan et al. (1996, 1997a) model of quiescent BH SXTs is the low radiative 
efficiency of the ADAF. In these flows, the bulk of the viscously dissipated energy is stored in the gas and 
advected with the flow into the black hole (Narayan & Yi 1995b; Abramowicz et al. 1995; Narayan et 
al. 1996, 1997a). This explains the unusually low luminosities of BH SXTs. By contrast, in NS SXTs 
all the advected energy is expected to be radiated from the neutron star surface, resulting in a much 
higher radiative efficiency of the accretion flow even in the presence of an ADAF (Narayan & Yi 1995b). 
Motivated by this fundamental distinction between black hole and neutron star systems, Narayan, Garcia & 
McClintock (1997) and Garcia et al. (1998) compared the outburst amplitudes of BH SXTs and NS SXTs 
as a function of their maximum luminosities and showed that the observations reveal systematically lower 
relative luminosities in BH SXTs. They argued that this constitutes a confirmation of the presence of an 
event horizon in BH SXTs. The argument was challenged by Chen et al. (1998). 

In this paper, we attempt to develop a physical understanding of the difference in quiescent luminosities 
between NS and BH SXTs. We first show, in §2, that there is indeed a significant difference in the observed 
quiescent luminosities of the two classes of objects, in contrast to the claim of Chen et al. (1998). We then 
use binary evolution models in §3 to estimate mass transfer rates in SXTs. We combine these estimates 
with the ADAF+thin disk accretion scenario, taking into account the presence of an event horizon in BH 
systems and a reradiating surface in NS systems, to determine the expected X-ray luminosities in quiescence 
(§4). The model predictions agree well with the observations of quiescent BH SXTs. However, the model 
substantially overestimates the luminosities of quiescent NS SXTs. In §5 we show that the NS SXT data 
can be reconciled with the predictions of the model if we take into account the "propeller effect" (Illarionov 
& Sunyaev 1975), whereby the magnetosphere of a rapidly rotating neutron star prevents much of the 
accreting material from reaching the surface of the neutron star. In §6 we show that an ADAF model with 
a moderate wind and a somewhat less efficient propeller (in NS SXTs) , is also consistent with the observed 
quiescent luminosities of BH and NS SXTs. Finally, in §7 we discuss possible limitations and extensions of 
this work and in §8 we summarize the main results. 

2. Observations 

In Table [|, we list key parameters of several NS and BH SXTs: the orbital period P rb, the distance 
D, the quiescent X-ray luminosity L m - ln in the 0.5 — 10 keV X-ray band, and the mass mi (in solar units) of 
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the compact primary. Since our main interest is in the quiescent emission, we list only systems for which 
there exist reliable measurements of L m ; n . The values of .L mm are taken from Narayan et al. (1997b) and 
Garcia et al. (1998), except for the BH SXT H1705-250 (see below). 

In selecting the systems listed in Table [j], we were careful to avoid any possible confusion between BH 
and NS SXTs, since this would make a comparison between the two sets of objects less reliable. Thus, 
we limit our sample to the eight BH candidates that have firm dynamical lower limits on the mass of the 
primary (e.g. McClintock 1998). (Note that this is not true for several so-called BH "candidates" of Chen 
et al. 1998, which seriously weakens their arguments.) Similarly, there is firm evidence that the five systems 
listed in Table | as NS SXTs contain NS primaries, based on the detection of type I X-ray bursts (e.g. 
Narayan et al. 1997b; Chen et al. 1997). 

Here we report a new and improved luminosity limit for H1705-250 (=Nova Oph 77) based on the ASCA 
observation of 18 September 1996, which we extracted from the HEASARC archive. The exposure time was 
31 ks, but unfortunately the target was near the chip boundaries which limits the effective area. The 4cr 
upper limit is 0.9 x 10~ 3 c s -1 , which corresponds to an X-ray flux (1 — 40 keV) < 1 x 10~ 13 ergs s cm~ 2 . 
Assuming a distance of 8.6 kpc, L x < 0.9 x 10 33 ergs s _1 . 

Luminosities can be expressed in units of the Eddington luminosity, 

L Edd = 1.25 x 10 38 mi erg s" 1 , (1) 

where mi, the mass of the compact object in solar units, is listed in Table |l|. For a standard radiative 
efficiency of 10 %, the Eddington luminosity corresponds to a mass accretion rate of 

M Edd = 1.39 x 10 18 mi g s -1 . (2) 

In this section, we use both absolute luminosities and Eddington-scaled luminosities^ since it is not 
clear which is the more appropriate quantity for comparisons. Our calculations in §3 suggest that the mass 
transfer rates in NS and BH SXTs are similar when expressed in units of MEdd) especially at short orbital 
periods. Thus, Eddington-scaled luminosities might be more appropriate to compare the two classes of 
objects at short periods (where most of the data lie). However, this argument is not very strong. Therefore, 
we also show absolute luminosities, expressed in units of 10 38 erg s" 1 (which is roughly the luminosity 
expected from an object accreting at 10 18 g s _1 with a radiative efficiency of 10 %). Finally, the ratio 
^mm/imax (which is independent of the distance D) is used for comparison with previous investigations of 
the luminosity difference between BH and NS SXTs. The values of L max , the outburst peak luminosity, are 
taken from Garcia et al. (1998). 

In Figure [l], we show the quiescent luminosities of the BH and NS systems listed in Table |l| as a function 
of their orbital periods P or b- The open circles correspond to NS SXTs and the filled dots correspond to BH 
systems. Luminosity upper limits (indicated by downward arrows) are shown for five SXTs, all of which 
are BH systems. The orbital period of the NS SXT H1608-52 is uncertain (see Table [l]) and its location in 
Figure [j] is indicated by dashed circles. Among the undetected SXTs, we choose to include in our sample 
only those systems that have been observed for more than 10 ks (see Table [j] herein, and also Table 1 in 



6 Note that Eddington-scaled luminosities are somewhat uncertain if the primary mass mi in a system is not well known 
(Eq. Jl]]). This is for instance the case of the BH SXT J1655-40. Orosz &z Bailyn (1997) estimate that mi = 7 in this system, 
while Phillips et al. (1999) claim that mi could be as low as 4.1 
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Narayan et al. 1997b), which therefore have flux limits < 3 x 10 13 ergs s 1 with current X-ray satellites 
(0.5-10 keV; ROSAT/ASCA). 

The argument for using P or b along the horizontal axis of Figure |l| is as follows. For any binary system 
with a low mass secondary, and transferring mass via Roche-lobe overflow, the density of the secondary 
essentially determines P or b- At a given P or b, a BH SXT and a NS SXT will have similar secondaries, so 
that the mass transfer characteristics are likely to be similar. Thus a reliable comparison of the radiative 
efficiencies of their accretion flows would be possible. In contrast, if we were to compare the quiescent 
luminosities of a BH SXT and a NS SXT with quite different P or b, a difference in the mass transfer rates 
(see §3) could mask actual differences in the radiative efficiencies of the accretion flows. This point, which 
is the motivation for the calculations in §3, was first emphasized by Lasota & Hameury (1998). 

Figure |l| strongly suggests that BH SXTs are fainter than NS SXTs. This is especially true when 
the data are plotted in Eddington units (Fig. ^a), but it is reasonably convincing even in absolute units 
(Fig. |l|b) or in terms of the ratio L m i n / L max (Fig. |l]c). Note that all five of the NS SXTs have been detected 
in quiescence and have well measured values of L m [ n , while the majority of BH SXTs have only upper limits 
on Z/min- In our view, this is significant and underlines the effect we are claiming. 

To make a clearer comparison between BH and NS SXTs with the existing data, we isolate in Figure ^ 
those systems for which L m j n and P or b are well determined (see also the very similar Fig. 3 of Lasota & 
Hameury 1998). Although there are only a few systems in this figure, we feel that we can conclude that (at 
least in this small sample) BH SXTs in quiescence are less luminous than NS SXTs. The difference is as 
much as two orders of magnitude (in Eddington units) for short period systems. 

It is important to stress here the excellent upper limit on the quiescent luminosity of the BH SXT 
J0422+32. This upper limit is nearly equal to the very low measured quiescent luminosity of A0620-00. 
Thus, there are now two BH SXTs that are two orders of magnitude less luminous, in Eddington units, 
than quiescent NS SXTs with comparable P or b- It is important to carry out more sensitive observations of 
those BH SXTs which have only upper limits currently. 

Narayan et al. (1997b) explained the difference in quiescent luminosities of NS SXTs and BH SXTs 
as due to the presence of a hard surface in NS SXTs vs an event horizon in BH SXTs. The presence of a 
surface causes the radiative efficiency in NS SXTs to be high, while in BH SXTs most of the dissipated 
energy is advected into the event horizon of the BH, without producing observable emission. Narayan et al. 
argued that this constitutes observational evidence for event horizons in black holes. In the remainder of 
the paper, we expand on this idea and construct a physical model of accretion in NS and BH systems to 
explain quantitatively the observations summarized in Figure ^. In §3 we estimate the mass transfer rate in 
SXTs, and in §4 and §5 we describe an accretion model that is able to reproduce the observed emission in 
the 0.5 - 10 keV band. 



3. Mass Transfer Rates Predicted By Binary Evolution Models 

The theory of binary evolution relies on the Roche-lobe model for the description of mass transfer 
(Frank, King & Raine 1992). Depending on which mechanism drives the mass transfer, whether it is loss of 
orbital angular momentum through gravitational radiation and magnetic braking, or expansion of the donor 
as it evolves away from the main-sequence, the binaries are classified as j-driven or n-driven systems. There 
is a 'bifurcation' orbital period Pbif separating the two classes such that for P or b > Pbif we find n-driven 
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systems whose orbital periods increase with time, and for P orb < Pbif we find j-driven systems whose P orb 
decrease with time. Although there is no unique value of Pbif, because it depends on the donor mass and 
the strength of angular momentum losses, estimates range from 0.5 — 2 days (e.g. Pylyser & Savonije 1988; 
King, Kolb & Burderi 1996). 

The exact form of the magnetic braking law is not well established (see e.g. Kalogera, Kolb & King 
1998). Moreover, Menou, Narayan & Lasota (1998) showed that there is strong circumstantial evidence 
that magnetic braking (MB) is weak in BH SXTs. In the following, for simplicity, we neglect the influence 
of MB in BH and NS SXTs. We will see in §4.2 that any contribution to the mass transfer by MB would 
only strengthen the argument that NS SXTs shed mass rather than accrete it. 

Following King et al. (1996), the mass transfer rates in j-driven systems with gravitational radiation 
(GR) and n-driven system with secondary expansion (EXP) is given by: 

M T « 1.27 x 10 14 mim|(mi + m 2 )- 1/3 (P orb /ld)- 8/3 g s" 1 , P orb < P bif (GR), (3) 
« 2.54 x 10 16 mi- 47 (P orb /ld) a93 g s"\ P orb > P bif (EXP), (4) 

where to 2 is the mass of the secondary in solar units, P orb is scaled in units of 1 day. 

The secondary in j-driven binaries is usually a main-sequence star (see King et al. 1996 for possible 
complications) filling its Roche lobe. These stars obey the simple relation m 2 = 0.11(P orb /l hr) (e.g. Frank 
et al. 1992). The prediction for rriT in case of GR driven mass transfer then becomes 

M T w 8.85 x 10 M mi(mi + m 2 )~ 1 / 3 (P orb /ld)- 2/3 g s" 1 , P orb < P bif (GR). (5) 

2/3 

which scales with the mass of the primary as m x for small m 2 . Note this is not very different from 
Mt oc mi, thus motivating the use of the Eddington scaling of luminosity in Figs. |l| and || 

Equation (Q) for the case of secondary expansion is the linear limit of third-order polynomial fits 
(King 1988) to detailed evolutionary tracks of single stars computed by Webbink, Rappaport & Savonije 
(1983). Along an evolutionary track, the mass to 2 of the expanding secondary decreases as it transfers mass 
onto the primary. Since a variety of initial masses and orbital periods (at the onset of mass-transfer) are 
possible, ra% is essentially a free parameter in equation (^). However, a careful study of evolutionary tracks^ 
shows that the band of mass transfer rates predicted by equation (Q) for values of to 2 in the range 0.5 to 
1 reproduces the mass transfer rates in n-driven accreting binaries for the majority of their lifetime. This 
band therefore provides, in a "population synthesis" sense, the likely mass transfer rates one can expect in 
SXTs with orbital periods larger than the bifurcation period. 

We estimate the mass transfer rate expected at a given orbital period by simply adding the contributions 
from GR (Eq. ||) and secondary expansion (a band with m 2 = [0.5, 1] in Eq. 0). Figure |3| shows that this 
results in a relatively short, although plausible, bifurcation period ~ 10 hr. The estimates of My are shown 
separately for BH and NS SXTs, in units of 10 18 g s _1 . The mass transfer rates are typically similar, in 
absolute units, in BH and NS SXTs with long orbital periods. However, the mass transfer rate in BH SXTs 
is larger than in NS SXTs at short orbital periods (cf Eq. and its scaling with mi); the two are in fact 
quite close in Eddington units at these periods, where the majority of SXTs are found. Although this result 
neglects any contribution from magnetic braking, we think that it is significant enough to motivate the use 
of Eddington-scaled luminosities in Figs. E] and 0. Note that the comparable values of Eddington-scaled rhr 



7 See Verbunt & van den Heuvel (1995) for a semi-analytical model that allows one to compute tracks. 
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predicted for BH SXTs and NS SXTs at short periods contrast with the large Eddington-scaled luminosity 
difference of the two classes of objects (Figs. [| and ||). This in essence is the argument of Narayan et al. 
(1997b) and Garcia et al. (1998), and we confirm their claim. The sample is, however, small. 

Although the estimates shown in Fig. |^ rely on simple considerations that do not take into account 
the full complexity of mass transfer (in particular, substantial fluctuations of the mass transfer rates 
around the secular values used here may occur; e.g. Warner 1995), they probably provide better than an 
order-of-magnitude precision, which is sufficient for our discussion. Note, however, that this simple model 
is not likely to give an accurate value of Mt for the BH SXT J 1655-40 since this system seems to be in a 
special evolutionary state (crossing the Hertzsprung gap, e.g. Kolb 1998). 



4. Accretion Models 

Following Narayan et al. (1996, 1997a), we assume that the accretion flow is made of two components: 
an inner ADAF and an outer thin accretion disk. Since only the ADAF contributes to the emission of 
quiescent SXTs in the 0.5-10 keV band, we ignore the thin disk. Further, since the emission of the ADAF is 
primarily from regions close to the center, we do not concern ourselves with the precise value of the radius 
i?tr at which the transition between the ADAF and the thin disk occurs. 

Although the ADAF model was initially proposed for quiescent BH SXTs, there is no reason why it 
should not apply also to quiescent NS SXTs. In particular, if the transition of the accretion flow from a 
thin disk configuration to an ADAF configuration is triggered by the same local physics in the disk, e.g. a 
thermal instability in the upper layers of the disk (Shaviv & Wehrse 1986; Meyer & Meyer-Hofmeister 1994; 
Narayan & Yi 1995b), the process should occur independently of the nature of the central object. We note, 
however, that the presence of a reradiating surface in the case of a NS primary at the center of the ADAF 
may have significant effects on the ADAF. For instance, Narayan & Yi (1995b) showed that the critical 
accretion rate m cr jt in an ADAF is reduced in the case of a NS (compared to a BH) because the additional 
supply of soft photons (from the stellar surface) causes the ADAF to cool more efficiently. This is not an 
issue here since m in all the quiescent systems we study is well below m cr it. 

We assume conservative mass transfer. Consequently, if A^adaf is the rate at which mass is accreted 
during quiescence via the ADAF and Af accum is the rate at which mass is accumulated in the outer thin 
disk, the mass transfer rate satisfies Mt = A/adaf + -^accum- The proportion of mass accreted vs. mass 
accumulated in quiescence is, however, uncertain. Menou et al. (1998) estimated the relative importance 
of A/adaf and Af accum for most of the BH SXTs listed in Table [|. They estimated M accum by taking the 
integrated fluences of the outbursts of each source and assuming a radiative efficiency of ry = 0.1 (which 
provides the total mass accreted in the outburst), and estimated Madaf by fitting the quiescent emission 
with an ADAF model. Menou et al.'s work indicates that, in quiescent BH SXTs, roughly half the matter 
that is transferred from the secondary is accreted via the ADAF and the rest is accumulated in the thin 
disk, i.e. 

A/adaf -, Madaf g 1 

—. = q ~ 1, ■. = ~ -. (6) 

M accum y M T 1 + 5 2 

In the following models, we treat g as a free parameter, but we expect g to be of order unity both in BH 
SXTs and in NS SXTs. 

The version of the ADAF model used here is described in detail by Narayan et al. (1998a). It 
consistently includes adiabatic compressive heating of electrons in the energy equation (Nakamura et al. 
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1997) and it uses for the flow dynamics the general relativistic global solutions calculated by Popham & 
Gammie (1998). 

In previous studies involving ADAFs, most of the model parameters were kept fixed at the following 
values (see Narayan et al. 1998b for a discussion): o;adaf = 0.3 (the viscosity parameter in the ADAF), 
= P g as/-Ptotai = 0.5 (i.e. equipartition between the gas and the magnetic field), 7 = 13/9 = 1.44 (adiabatic 
index of the gas) and 5 = 10~ 3 (fraction of the viscous dissipation that goes directly into heating the 
electrons). We adopt these standard values of c^daf, P and 7 in the present models (see §6 for the effect 
of varying 8). We take Mt to be given by the predictions of binary evolution models (Fig. and obtain 
-Madaf via equation (^]) for a given assumed value of g. We calculate the spectral energy distributions of 
the models and numerically integrate the spectra to obtain the luminosities in the 0.5-10 keV band. The 
results are discussed in the following sections. 

4.1. ADAF models of BH SXTs 

Figure shows two bands of luminosities corresponding to ADAF models of BH SXTs. In the upper 
band Madaf = Mt, i.e. the entire Mt is accreted by the black hole via the ADAF, while in the lower band 
A/adaf = Mt/3 (or g = 1/2). The lower band is in satisfying agreement with the observed luminosities. 
The fact that we obtain agreement with a reasonable value of g means that the modeling of the quiescent 
X-ray emission of BH SXTs with the ADAF model is consistent with the predictions of binary evolution 
models. The increase in the luminosity with increasing P or b is also in good agreement. One should stress 
again, however, that the agreement for the BH SXT J1655-40 is unexpected if the system is crossing the 
Hertzsprung gap (Kolb 1998). 

The bolometric radiative efficiency in the ADAF models constructed here is typically ~ 10~ 3 compared 
to ~ 10 _1 for standard thin disk accretion. The X-ray radiative efficiency of these ADAF models is another 
two or three orders of magnitude lower since much of the ADAF luminosity comes out as synchrotron 
radiation in the optical. Note that these results are sensitive to various details of the ADAF models, and in 
particular to the choice of «adaf and f3. For instance, we find that, if «adaf = 0.1, we need g ~ 1/4 (i.e. 
a fraction ~ 1/5 of the mass supplied by the secondary must be accreted via the ADAF) in order to obtain 
agreement with the observed quiescent luminosities of BH SXTs. The fraction goes down further to ~ 1/10 
if aADAF = 0.025. 

4.2. ADAF models of NS SXTs 

The presence of a hard surface on the accreting star is the key feature of ADAF models of NS SXTs. 
The large amount of energy advected in the ADAF, which is lost through the event horizon in BH SXTs, 
will be reradiated from the NS surface (Narayan & Yi 1995b). Narayan et al. (1997b) and Garcia et al. 
(1998) argued that this additional source of luminosity is the explanation for the systematically larger 
luminosity of quiescent NS SXTs as compared to quiescent BH SXTs (Figs [l] and ||). We now investigate 
this quantitatively. 

In the following, we assume that the thermal energy stored in the ADAF is reradiated as blackbody 
emission with a radiative efficiency i] ~ 0.1. The luminosity in the 0.5-10 keV band will be dominated by 
the emission coming from the NS surface, so that the additional emission from the ADAF can be ignored 
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(for simplicity). We also assume that the reradiation occurs from a small fraction / sur f of the NS. This is in 
agreement with observations of quiescent NS SXTs (e.g. Verbunt et al. 1994; Campana et al. 1998; but see 
Rutledge et al. 1998 and the discussion in §7). For now, we do not seek to justify the small emitting area 
of the NS. We simply assume for all our NS models that the area is 10 km 2 . However, our results are not 
seriously affected even if the radiating area is increased or decreased by, say, an order of magnitude (i.e as 
long as most of the blackbody emission is in the 0.5-10 keV X-ray band). The typical effective temperatures 
found are T off < 0.5 keV. 

The upper band in Fig. ^b shows the luminosities predicted by binary evolution models if g = 1/2, i.e 
a third of the mass transferred by the secondary is accreted onto the NS (as in BH SXTs). We see that 
the luminosities predicted are much larger than those observed, by ~ 3 orders of magnitude (except for the 
system EXO 0748-676; see §7 for a discussion of this system). Since any contribution to the mass transfer 
by magnetic braking would only increase the mass transfer rates expected (in short orbital period systems), 
the discrepancy between the observed and predicted luminosities would be even more serious if magnetic 
braking were effective. 

One possibility for the discrepancy is that most of the quiescent flux from NS SXTs is emitted outside 
the 0.5-10 keV X-ray band. This explanation seems, however, rather unlikely since this emission would 
be seen somewhere else in the spectrum, say in hard X-rays or soft 7-rays. Aql X-l has been observed 
in quiescence in hard X-rays and it is clear that the energy does not come out in the 10-100 keV band 
(Campana et al. 1998). In addition, although some theoretical models of boundary layers between an 
accretion flow and a NS predict substantial deviations from blackbody emission at higher energies (e.g. 
Shapiro & Salpeter 1975, Turolla et al. 1994), there are no models predicting that such a large fraction 
(~ 0.999) of the emission would come out in hard X-rays or 7-rays. 

We must therefore consider the possibility that only part of the mass supplied by the secondary 
actually reaches the surface of the NS. The lower band in Fig. ||b shows that the luminosities of quiescent 
NS SXTs can be accounted for (except EXO 0748-676) if only ~ 0.1% of the transferred mass reaches the 
NS surface. This is of course a surprisingly small fraction.^. In the following section, we argue that the 
propeller effect can explain the fraction if combined with accretion via an ADAF, while in §6 we consider 
the possibility that winds from ADAFs could further reduce the amount of mass reaching the NS surface. 



5. Propeller Effect in Quiescent Neutron Star SXTs 

The propeller effect in accreting neutron stars was initially proposed by Illarionov & Sunyaev (1975) 
to explain the existence of long-period X-ray pulsars, and it was later developed in more detail by Davies 
and Pringle (1981; see also Wang & Robertson 1985; Stella et al. 1986). Although the details vary, the 
basic mechanism of the effect relies on the presence of a centrifugal barrier at the rotating magnetosphcrc 
of a rapidly spinning NS. 

Inside the magnetosphere, the accreting gas is forced to follow the magnetic field lines of the NS 
since, by definition, in this region magnetic forces dominate the flow dynamics. The fate of the gas is then 
determined by the relative magnitudes of the magnetospheric radius, i? m , and the corotation radius, R co , 



8 Although this value might simply be interpreted, in a scenario without ADAF, as the accretion rate at the inner edge of 
an unsteady quiescent disk (M oc R 2,5 ~ 3 in such a disk; e.g. Cannizzo 1993), we do not favor this interpretation because the 
propeller effect should operate in this case and prevent any mass from reaching the NS surface (see §5) 
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defined by = (GM/i?^) 1 ' 2 , where SI* is the angular velocity of the NS (and the magnetosphere) . If 
i? co < Run the magnetosphere rotates so fast that the centrifugal force at R m is larger than the force of 
gravity and hardly any of the gas can be accreted onto the NS. This is the propeller effect. The efficiency of 
the propeller effect is not well understood (see Davies, Fabian, & Pringle 1979, and Davies & Pringlc 1981). 
It is also not clear whether the accreting matter is ejected from the system or is merely accumulated in a 
boundary layer around the magnetosphere (Wang & Robertson 1985). However, from the simple physical 
argument outlined above, it is clear that normal accretion onto the neutron star is inhibited during the 
propeller phase. 

The first direct observational evidence for the existence of the propeller effect in accreting neutron stars 
was reported by Cui (1997) in two X-ray pulsars, GX 1+4 and GRO J1744-28. Cui showed that in both 
systems X-ray pulsations ceased during a period of low X-ray flux, which he interpreted as the result of a 
decrease in the mass accretion rate and a corresponding increase in the magnetospheric radius to beyond 
the corotation radius. Campana et al. (1998) and Zhang et al. (1998) argued that the propeller effect was 
also seen in a NS SXT, Aql X-l, during its most recent outburst; the system showed an abrupt decrease 
of the X-ray flux, accompanied by considerable hardening of the X-ray spectrum. Both sets of authors 
proposed that the hard spectrum originates just outside the magnetosphere, where the gas becomes very 
hot due to the action of the propeller (see also Cui et al. 1998). 

Observations of Aql X-l show that even in quiescence the blackbody spectral component does not 
disappear entirely implying that a small amount of material continues to accrete onto the star, contrary to 
simple models of the propeller (Verbunt et al. 1994; Stella et al. 1994). To resolve this problem, Zhang et 
al. (1998) proposed that accretion in quiescent NS SXTs occurs via a quasi-spherical ADAF, rather than a 
thin disk. This allows some material to accrete near the poles, thereby bypassing the centrifugal barrier. In 
this section we develop this idea quantitatively. 

5.1. Accretion Geometry in Quiescence 

In the following, we assume for simplicity that the NS spin axis is perpendicular to the binary orbital 
plane and that the NS spin and magnetic axes are aligned. 

5.1.1. Magnetospheric Radius 

In the presence of a spherical accretion flow (e.g. Bondi 1952), the magnetospheric (or Alfvcn) radius, 
i? m , is usually defined as the radius at which the magnetic pressure, P mag = B 2 /8ir, due to the neutron 
star magnetic field balances the ram pressure, P ram = P v r/2, of the accreting gas, where v r is the radial 
infall velocity of the gas (e.g. Frank et al. 1992). The contribution from the thermal pressure of the gas 
is generally neglected when considering Bondi-type flows which are highly supersonic (free- falling) . An 
ADAF, on the other hand, is subsonic down to a few Schwarzschild radii (e.g. Narayan, Kato, & Honma 
1997; Chen, Abramowicz & Lasota 1997; Abramowicz, Chen, Granath & Lasota 1996) and the thermal 
pressure term must be taken into account. 

The polar structure of ADAFs was described by Narayan & Yi (1995a) in the self-similar approximation. 
They showed that in the limit of strong advection (/ a dv ~ 1), the gas density and temperature (and 
therefore the thermal pressure of the gas) are essentially independent of the polar angle 0; i.e. these 
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quantities are nearly constant on spherical shells centered on the accreting object. On the other hand, 
the radial infall velocity is smaller than the thermal sound speed in the equatorial plane and decreases 
towards the poles. Clearly then, the thermal pressure of the accreting material is the dominant term in 
calculating the magnetospheric radius. Since this pressure is constant with 9, the magnetosphere will be 
nearly spherical. 

Despite the difference between a Bondi spherical accretion flow and an ADAF, the value of R m is 
roughly the same in the two cases. In the former, essentially all the gravitational energy released during 
accretion goes into bulk kinetic energy of the flow, while in ADAFs, the energy is stored as thermal energy 
of the gas. Thus, the ram pressure in Bondi flows must be of the same order as the thermal pressure in 
ADAFs. The magnetospheric radius in quiescence then takes the same form (Frank et al. 1992): 

/ ' \ ~ 2/? 

R mq = 6.45 x 10 5 -r— m- 3/7 B 8 4/7 4 2 s /7 6 cm - (?) 

\A/Edd J 

where B s is the NS surface magnetic field in units of 10 8 G and i?NS,6 is t ne NS radius in units of 10 6 cm. 

For the remainder of the paper, R refers to the radius in physical units and r refers to the radius in 
Schwarzschild units. 



5.1.2. Accretion in the Propeller Phase 

The propeller effect in binary systems is generally discussed in the context of thin disk accretion. In 
this model, once the propeller becomes effective, no matter is able to reach the NS surface because the 
centrifugal acceleration acts equally against gravity for all the matter located in the disk orbital plane. 

In a spherical flow, however, the centrifugal acceleration at the magnetospheric radius acting on a 
parcel of gas, accreting at a polar angle 9 from the spin axis, is equal to A c = fi^i? mq sin 9, where is the 
angular speed of the magnetic field lines anchored in the NS. The direction of this force is perpendicular 
to the spin axis of the NS (and parallel to the orbital plane), so that the component of A c along the 
radial direction is simply A c sin 9 = f2ji? mq sin 2 9. The maximum polar angle below which the gravitational 
acceleration (= i£mq^if(-Rmq)) wins over the centrifugal force is given by 

sin l?o = ^ — — . (8) 



For 9 <C 1 this simplifies to 



9 * ^g^- (9) 



Note that the residual force acting on all parcels of gas, even those which do not overcome the 
centrifugal barrier, tends to direct them toward the orbital plane. In this simple formulation, we do not 
consider possible complications due to field line orientation effects. 
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5.1.3. Fraction of Mass Reaching the NS Surface 

Only matter accreting between 6 = and 9 = 9 can overcome the centrifugal barrier and reach the 
surface of the NS. The mass accretion rate onto the star is then 

M NS = -2 / 2nR sin 9pv r Rd6, (10) 



o 



where p = p{R, 9) ~ p(R) is the gas density and v r = v r (R, 9) ~ v r (R) sin 2 9 is the radial infall velocity at 
angle 9 (positive outward), as given by the self-similar solution of Narayan & Yi (1995a). The factor of 2 in 
equation ([To]) is introduced to account for the two polar caps. The fraction of the total mass accretion rate 
Madaf in the ADAF that reaches the NS surface is then 



A%s 2 Ia° 2lrR2 sin 0p(R)v r (R) sin 2 9d9 2 f°° sin 3 



Madaf 2 f* /2 2vri? 2 sin 9p(R)v r (R) sin 2 9d9 2 f* /2 sin 3 I 

(11) 



3 ^4 3 f ^x(i?mq) X 



I A' 



where the last two steps have been calculated in the limit of small 9q. For a neutron star with a spin 
period of 2 ms, a magnetic field strength of 3 x 10 s G and accreting at Madaf = 10 _3 MEdd, we find 
/ acc ~ 0.75 x 10 -3 . Note that if the radial velocity does not vary as sin 2 6 but is more nearly constant with 
9 (as in Bondi accretion), then / acc would be larger (e.g., if v r (R,9) is independent of 9, / acc oc 9q only). 
On the other hand, if the accretion flow has a toroidal morphology with empty funnels along the rotation 



axis, / acc would be much lower than the expression given in equation (11), perhaps even zero. There is thus 
considerable uncertainty in the value of / acc . 



5.1.4- Fraction of the NS Surface Emitting Radiation 

We assume that the gas accreted during the propeller phase follows the magnetic field lines down to 
the NS polar caps, where its kinetic and thermal energy is converted to radiation. Dipolar field lines satisfy 
the parametric equation (e.g. Frank et al. 1992): 

i? = Csin 2 6>, (12) 

where C is constant for a given field line. Thus, the field lines intersecting the magnetospheric radius at 
9 = 9q, emerge from the NS surface at a polar angle 9s given by: 

sin 2 9 S = R NS ^ 

Sin 2 8 Rmq 

The accreting material reaches the NS surface at 9 < 8s, and therefore the fraction of the neutron star 
surface that re-radiates the accreted energy is 

_2j^7^n9d9_ 9J 8jR^s 

/surf - /2 _ — _ — _ „ , 



2 r;/ 2 27 rsin^ - ? R 



mq 



where we have made use of the fact that 6s < Go <C 1- Note the direct dependence of / sur f on the NS 
magnetic field strength B through R mq . For a neutron star with a spin period of 2 ms, a magnetic field 
strength of 3 x 10 8 G and accreting at Madaf = 10~ 3 MEdd, we find / sur f — 3 x 10~ 3 . 
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5.2. Equilibrium Spin Frequency 

Equations ( pi] ) and (|l4|) demonstrate that in quiescence the mass accretion rate onto a NS and the 
emitting area, which determine the total luminosity and spectrum of the system, depend strongly on the 
NS angular rotation speed f2* and the field strength. Although Q± is not known for most NS SXTs, we 
discuss in this section how it can be estimated using our knowledge of the accretion history of an SXT and 
of the interaction between the magnetosphere of the NS and the accretion flow. 

During the propeller phase, in quiescence, the NS ejects most of the accreting gas that reaches its 
magnetosphere. The ejected material leaves the system with higher specific angular momentum than it 
had coming in, which results in an effective spin-down of the NS. However, NS SXTs experience outbursts 
during which the mass accretion rate onto the NS is orders of magnitude higher than in quiescence. High 
M causes the magnetosphere to shrink inside the corotation radius so that all the accreting gas is able to 
reach the NS surface. The mass and angular momentum of the accreting gas is added to the NS, which 
causes it to spin- up. A NS SXT that accretes for a fairly long time is expected to reach an equilibrium spin 
frequency 0* jeq such that the spin up during outburst is balanced by the spin down during quiescence. 

At the onset of mass transfer, the spin period of the NS in an SXT is likely to be different from the 
equilibrium value. In standard Ghosh & Lamb (1979) type models (for steady disks), the timescale r oq on 
which a typical NS reaches the equilibrium spin period is usually small (> 10 5 yr) compared to the lifetime 
of the LMXB (10 s — 10 9 yr; see, e.g., Henrichs 1983 for a review). Although T eq may be a little larger for 
transient accretion, it is reasonable to assume that most of the neutron stars in SXTs have rotation rates 
close to f2*, C q- It is straightforward to show that is basically unaffected by a single outburst or a single 
quiescent phase, because the time spent between two outbursts (~ a few years) is much shorter than r cq . 
Therefore, it is not necessary to consider the small jitter in f2* during the quiescence-outburst cycle. In 
the following, we simply assume that once a NS SXT has reached spin equilibrium, then = f2* jCq both 
during outburst and quiescence. 



5.2.1. Idealized Model for the Spin-up and Spin-down Torques 

To determine the equilibrium spin frequency of a NS in a transient binary system we need a theoretical 
description of the interaction between the NS and the accretion flow. We begin with the outburst phase 
during which accretion occurs via a thin disk and the torque exerted by the accreting gas on the NS is 
relatively well understood. 

Despite some modifications, current models describing the interaction between a NS and a magnetically 
threaded thin accretion disk remain essentially identical to the model initially proposed by Ghosh & 
Lamb (1979). In this model, the spin-up torque on the NS results from two contributions. First, the 
gas that follows magnetic field lines and reaches the NS surface gives its angular momentum to the NS 
(spin-up). Second, the interaction between the magnetic field lines and the threaded thin disk beyond the 
magnetospheric radius results in a positive torque on the NS for radii where the field lines rotate slower 
than the local Keplerian angular speed of the gas. The lines that thread the disk further out give a (smaller) 
negative torque on the NS since these field lines rotate faster than the local Keplerian angular speed of 
the gas (Gosh & Lamb 1979). Recent detailed numerical simulations by Daumerie (1996) show that the 
overall contribution resulting from the interaction of the accretion flow and the magnetic field lines beyond 
the magnetospheric radius is nearly equal to the contribution to the spin-up from the gas reaching the 
NS surface. Further, Daumerie finds that the torque varies roughly linearly with the fastness parameter 
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w = si*/si K (Rm). 

Using this work as a guide, we propose the following idealized formula for the torque: 

j = 2MR 2 m fl K (R m ) 

where the factor of 2 appears because of the two nearly equal contributions to the spin-up mentioned above. 
Equation ( |l5|) assumes that the fastness parameter at equilibrium, uj c1 -h, is exactly unity. This is consistent 
with the recent arguments of Wang (1995; see also Wang 1987). 

The torque in the propeller regime, when f2* > £li<-(i? m ), is not as well constrained. By analogy with 
Ghosh & Lamb type models, we assume that the torque arises from two nearly equal contributions. One 
contribution is the negative torque applied to the NS when the propeller expels gas at the magnetosphere 
and the other is the negative torquef] due to the interaction of the magnetic field lines with the accretion 
flow beyond the magnetospheric radius. If the gas expelled leaves the magnetosphere with an angular 
momentum corresponding to the angular rotation speed of the NS, viz. (which corresponds to an 
efficient propeller), then the negative torque on the NS in the propeller phase is 

j ~ -2Mi^Q*, (16) 

where M is the rate at which mass reaches the magnetosphere. (The small positive contribution from 
the very small fraction of mass accreted on to the NS is neglected.) This expression is in agreement with 
equation (^|) in the limit of f2* ^> Qk {Rm)- Therefore, we can use the prescription given by equation ( |l5| ) 
to describe the interaction between the NS and the accretion flow throughout the outburst-quiescence cycle. 



1 - 



^k(Rtd) 



(15) 



5.2.2. Spin Equilibrium 



We define AJ a — J At to be the integrated spin-up torque over an outburst of duration At Q , 
and AJ q = J q At q to be the integrated spin-down torque during a propeller phase of duration At q . 
The recurrence time of the SXT is At + At q . We assume, for simplicity, that the two phases can be 
approximated as bimodal, with the outburst phase occurring at a typical accretion rate M and the 
quiescent phase having a typical lower accretion rate M q . These accretion rates define, in turn, the values of 
the magnetospheric radius in outburst,^] i? mo , and in quiescence, i? mq , used in equation (|l^) to determine 
J and J q . Since the accretion rate from the peak of the outburst is known to decrease gradually with time, 
the value of M Q used to determine R mo should be interpreted as a mean over the outburst. 

By definition, the NS angular rotation speed at equilibrium f2* )e q satisfies AJ a = AJ q . Using 
equation ([L5|), we find 



ft*,, 



Q,K{Rn 



Rn 



Rn 



3/2 



1 



1/2 



(17) 



9 This contribution to the torque is always negative in the propeller regime because SI* > Qpc(R) for all R > Rmq- 

10 Following Frank et al. (1992), we use a value for the magnetospheric radius of the thin disk that is half that used for 
spherical accretion (Eq. M) 



-14- 



where g is the ratio M q At q /M At , i.e. the total mass reaching the magnetosphere during the ADAF- 
propeller phase divided by the total mass accreted during outburst. This definition of g is equivalent to the 
definition of g given in equation (|J) if all the mass accumulated in the disk during quiescence is accreted 
onto the neutron star during outburst.^ For the remainder of the paper, we assume g — g. Note that the 
spin equilibrium period defined by equation ( |i~7| ) differs from the spin equilibrium period defined for steady 
accretion onto a NS and usually discussed in the literature (e.g. Henrichs 1983). 



5.3. Propeller Regime at Spin Equilibrium 

Equation ( |l7j ) allows us to determine the equilibrium spin period P sp i n as a function of the parameter 
g = A/adaf / -^accum ■ This, in turn, allows us to estimate the quantities / acc (Eq. |ll|) and / sur f (Eq. Jd]) 
which are relevant for determining the observational properties of quiescent NS SXTs. In the following 
calculations, we assume for definiteness that the NS accretes at the Eddington rate (M = MEdd) in 
outburst and that the mass transfer accretion rate in quiescence is Mr = 10~ 2 ' 5 AiEdd (cf Fig. ||). However, 
our conclusions do not depend crucially on these assumptions. 

Figures |^ and ^| show the variations of M^/Mt = /acc x 5/(1 + <?)) /surf and -P S pm with g. As Madaf 
decreases, the magnetosphere becomes more extended and the propeller effect becomes more efficient. The 
solid lines show the results for a NS with a magnetic field strength B = 10 s G, while the dashed and dotted 
lines show the cases B = 10 9 and B = 10 10 G, respectively. Recent observations suggest that accreting 
neutron stars in low mass binaries have rather low magnetic fields, typically B < 10 9 G (e.g. White & 
Zhang 1997). 

Figure || can be related to Fig. ||b, where we found that M^s/Mt ~ 10~ 3 was needed in order to 
explain the very low observed luminosities of quiescent NS SXTs. Such a value of M-^s/Mt is predicted by 
the propeller model if g ~ 1/5. This value of g is reasonably close to the value g ~ 1/2 that we obtained for 
BH SXTs. Note that the dependence of M-^s/Mt on the magnetic field B should cancel out through the 
ratios Rm q /R mo in equation (|l^). For low magnetic field strengths, however, the inner edge of the disk in 
outburst is not given by half the value in equation (R) but is fixed at the last stable orbit at 3 Schwarzschild 
radii (for instance, if M a > O.lMEdd and B < 10 9 G). Consequently, there is a residual dependence of 
M NS /M T on B for small B. 

Figure |6^i shows that the fraction / sur f of the NS surface that emits in quiescence depends much more 
strongly on B than / acc does (Fig. ||), as expected from equation Jl4|). For values of g ~ 1/5 and low 
magnetic fields (B < 10 9 G), the emitting surface is typically between 1 and 10 km 2 . Observations of 
quiescent NS SXTs suggest an emitting surface ~ 1 km 2 (e.g. Verbunt et al. 1994; Campana et al. 1998; 
but see §7). Given that the surface could very well be underestimated by the standard techniques (e.g. 
Lewin, Van Paradijs & Taam 1993; Rajagopal & Romani 1996; Zavlin et al. 1996; Rutledge et al. 1998), 
values of ~ 1 — 10 km 2 appear in reasonable agreement with the observations. We note, however, that this 
result depends quite crucially on the assumed alignment of the NS spin and magnetic axes. This is a major 
uncertainty in the calculations. 

Figure ^b shows that P sp i n is even more sensitive to the value of B. For weak magnetic fields (B < 10 9 



1 A possible situation where the mass accreted onto the neutron star during outburst is less than the mass accumulated in the 
disk during quiescence is if the depletion rate of the disk during outburst is super-Eddington while the accretion rate onto the 
neutron star is Eddington-limited. Even in that case, Eq. (II 7|) remains valid because 0* >eq depends on the ratio M q At q /M Af Q 
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G), the values of the equilibrium spin period are in the range 1 — 10 ms, which is in agreement with values 
recently discussed in the literature for accreting neutron stars in low mass binaries (e.g. White & Zhang 

1997) . 

There are several sources of uncertainty in our simplified model for the propeller, which are discussed 
in §6.2 and §7 below. The proximity of the value g ~ 1/2 required to explain the quiescent luminosities of 
BH SXTs to the value g ~ 1/5 required to explain the luminosities of NS SXTs (with a propeller) suggests, 
however, that roughly the same geometry of accretion (outer thin disk and inner ADAF sharing the mass 
supplied by the secondary roughly equally) can account for the properties of both classes of objects. The 
presence of a hard surface and a magnetic field associated with the spinning neutron star are two major 
ingredients of the model. The NS SXTs Aql X-l (B ~ 10 8 G, F spin ~ 2 - 3 ms, emitting area ~ 1 km 2 ; 
White & Zhang 1997; Zhang et al. 1998; Campana et al. 1998) and SAX J1808.4-3658 (B < 2 x 10 s G, 
P spin = 2.49 ms; Wijnands & van der Klis 1998; Chakrabarty & Morgan 1998, hereafter CM98) are two 
interesting candidates that could allow further testing of the predictions of our ADAF-propeller scenario. 

We note that CM98 inferred a mass transfer rate in SAX J1808.4-3658 « 1CT 11 M Q yr -1 , which is 
roughly a factor ten less than the value for a 2 hr orbital period NS SXT according to our Fig. ||. These 
authors find that the estimate is consistent with mass transfer being driven by angular momentum losses due 
to gravitational radiation, if the secondary star has a mass ni2 = 0.05. The values shown in Fig. [|, however, 
were deduced assuming that the donor star is a main sequence star which satisfies to 2 = 0.11(P or b/l hr) 
(§3). This is apparently not the case in SAX J1808. 4-3658, perhaps because of the fact that the companion 
star is strongly irradiated (CM98). 

6. Winds from ADAFs and the Efficiency of the Propeller Effect 

The gas in an ADAF could in principle escape to infinity via a wind because part of the flow may be 
unbound. Narayan & Yi (1994, 1995a) showed that the Bernoulli constant, a measure of the energy of the 
gas at infinity, is positive in their self-similar solution, and Narayan, Kato & Honma (1997) showed that this 
is also the case in a global solution. A simple way of taking this into account was proposed by Blandford & 
Begelman (1998; hereafter BB98) who derived self-similar ADAF solutions that include mass loss via winds. 

Our results in §4.1 suggest that reducing chadaf has only a small effect on the radiative efficiency 
of the accretion flow. Quataert & Narayan (1998) showed, however, that an efficient way of increasing 
the radiative efficiency of the ADAF is to include winds and simultaneously increase the value of the 
parameter 5, the fraction of the viscously dissipated energy that directly heats the electrons. They found 
that models with winds and high S are consistent with the current observations for the BH SXT V404 Cyg 
in quiescence and the Galactic Center source Sgr A*. Note that theoretical estimates of the value of 5 are 
highly uncertain. Quataert (1998) and Gruzinov (1998) argued that S could be small if the magnetic fields 
in the ADAF are fairly subthermal, but Bisnovatyi-Kogan & Lovelace(1997; see also Quataert & Gruzinov 

1998) argued that magnetic reconnection in the accretion flow could preferentially heat the electrons and 
lead to large values of 5 (~ 1). 

In §6.1 and §6.2, we investigate the following two questions. Can winds by themselves, i.e. without any 
propeller effect, explain the observed quiescent luminosities of NS SXTs? ff not, is a combination of winds 
and a propeller consistent with the observations? 
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6.1. ADAF Models of NS and BH SXTs with Winds 

Following BB98, we assume that the accretion rate in the ADAF scales in a self-similar way, i.e. 

m(r) = m(r tr ) (j-^j , (18) 

where m is the accretion rate scaled in Eddington units (m = M/M-^^), < p < 1, and we assume that 
the wind is effective over the entire radial extent of the ADAF out to r out = r tr , the transition radius. 

The efficiency of mass loss is not only determined by p, but also by the radial extent of the ADAF 
(Eq. Around a black hole, the ADAF extends from the transition radius (r tr ) down to the event 

horizon, and so the accretion rate onto the black hole is fn- m = m(rt r )(l/Vt r ) p . Around a neutron star, 
however, the ADAF has a minimum radius equal to the magnetospheric radius r m , so that the accretion 
rate at the magnetosphere is mi n = m(rt r )(r m /r tr ) p . If r m 3> 1, as is the case for r mq in quiescent NS 
SXTs, then for a given value of p, the presence of a magnetosphere around the NS significantly reduces the 
total mass loss relative to the BH case. The value of r m itself depends on the local accretion rate at the 
magnetosphere (Eq. Since the accretion rate at a given radius decreases with increasing p, a larger p 
implies both a larger r m and a smaller radial extent for the ADAF. 

We solve for the magnetospheric radius by combining equations (^) and (18) and find 

i 

r mq = [l.35 x m(r tr )~ 2 r 2 r p 5 4 i?^ 6 ] 7+2p , (19) 
where r mq is expressed in Schwarzschild units for a NS of 1.4 M Q (Eq Il9l is the generalization of Eq M for 



non-zero p). We have confirmed (by combining Eqs. |18| and |19|) that despite the reduction of the radial 



extent of the ADAF for larger values of p, increasing p leads to an effective reduction of m(r mq ) for any 
reasonable values of p, r tr and m(rt r ). 

For a standard NS with low magnetic field strength (B$ — 1, i?NS.6 — 1) and a typical accretion rate at 
the outer boundary of the ADAF m(r tr ) = 1/3 x 10~ 25 (§4.2; Fig. |||), we find r mq ~ 10 in the absence of a 
wind (p = 0). For a wind with p = 0.4, r tr = 10 4 ("intermediate" wind model, 5 = 0.3), we find r mq ~ 20 
and m(r mq )/m(rtr) — 8 x 10 -2 , while for a wind with p = 0.8, rt r = 10 4 ("strong" wind model, S = 0.75), 
we find r mq ~ 35 and rh(r mq ) / 'm(r tr ) ~ 1 x 10~ 2 . For these estimates, we have chosen r tr = 10 by analogy 
with the values usually inferred in quiescent BH SXTs (Narayan et al. 1996; Menou et al. 1998). Were r tr 
to be smaller, the overall mass loss in the wind for a given p would be smaller. Note that the ADAF cannot 
extend far beyond 10 4 Schwarzschild radii because, for a given m(r tr ), there is a maximum radius out to 
which the ADAF can exist (e.g. Esin et al. 1997; Menou et al. 1998). 

Following Quataert & Narayan (1998), we computed two examples of ADAF models assuming that 
the dominant effect of the wind on the structure and the emission properties of the ADAF is to reduce 
the density in the accretion flow; in addition, we assumed that the emission from the wind itself can be 
neglected relative to the emission from the accretion flow. 

Figure [7]a shows the quiescent X-ray luminosities (0.5-10 keV) of BH SXTs predicted by the 
ADAF+wind model for which rh(r tr ) is equal to 1/3 of the mass supplied by the secondary p = 0.4, 
rt r = 10 4 and S — 0.3. This model gives roughly the same X-ray luminosities as the previous no-wind 
ADAF model (§4.1). The actual mass accreting onto the black hole is reduced here, but nevertheless the 
luminosities are the same because of the increase in the value of S from 10~ 3 to 0.3 (see Quataert k, Narayan 
1998). 
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Figure [?]b shows the quiescent X-ray luminosities (0.5-10 keV) of NS SXTs predicted by the same 
ADAF+wind model, assuming that all the mass reaching the magnetosphere is accreted onto the NS surface 
(i.e. no propeller effect). Since the X-ray luminosity of the ADAF is still much less than the luminosity 
coming from the NS surface, the main effect of including the wind is to decrease the luminosity of quiescent 
NS SXTs in proportion to the reduction of m at the inner edge of the ADAF (Eq. fljsfl). By comparing 
figures and 0b, it is clear that the luminosity is smaller by roughly an order of magnitude, which is 
insufficient to explain the observations. Thus, the intermediate wind scenario must also invoke a propeller 
effect, though with a somewhat lower efficiency than that described in §5. Given the uncertainties in our 
propeller model (see §6.2), we see no objection to such a model. 

Only for the strongest wind models (p ~ 1, r tr = 10 4 ) does the luminosity difference between NS SXTs 
and BH SXTs predicted by the models reduce to the observed difference. For p — 1 and r tr = 10 4 (Bg = 1, 
-Rns,6 = l)j Eq- (HI) an d fll9| ) predict r mq ~ 58 and m^^^/m^^) ~ 6 x 10~ 3 , which is basically consistent 
with the small accretion rates required to explain the observed quiescent luminosities of NS SXTs (cf. 
Fig. [|b). But is this a reasonable model? Can the propeller effect be completely ineffective at preventing 
mass from reaching the NS surface? 



6.2. The Efficiency of the Propeller Effect 

If the neutron stars in NS SXTs have spin periods of typically a few ms, as suggested by recent 
observations (e.g. White & Zhang 1997; Zhang et al. 1998; Campana et al. 1998; Wijnands & van der Klis 
1998), their corotation radii are a few Schwarzschild radii. This is substantially smaller than the few tens of 
Schwarzschild radii inferred for r mq in §6.1. The propeller effect is therefore expected to act quite strongly 
in these systems (lllarionov & Sunyaev 1975). Note that the values of r mq quoted in §6.1, which correspond 
to a magnetic field strength B of 10 8 G, would be even larger for a larger B (eq. ]19]| ). Can the propeller in 
such a system be so ineffective that / acc ~ 1? 

We pointed out in §5.2 that there are large uncertainties in the values of / aC c, which come mainly from 
the strong dependence of f acc on the NS rotation speed and the magnetospheric radius: / acc oc i7~ 4 i?~ q 
(the dependence on i? mq also affects the size of the NS emitting area, since / sur f oc R^ q )- The latter source 
of uncertainty is perhaps more important, since while can be determined by direct observations, R mq 
must be deduced from theoretical arguments. 

Although the value of R mq used in our work (Eq. [Q) is standard, it is based on simple dimensional 
arguments. Given the strong dependence of / acc on i? mq , an error, for instance, of only a factor of 1.5 in i? mq 
translates to an error of more than an order of magnitude in / acc . Recent work by Psaltis & Chakrabarty 
(1998) shows that the standard scaling used for R mq (for disk-magnetosphere interactions) is not consistent 
with observations in at least one weakly magnetized NS SXT (SAX J1808. 4-3658). Consequently, we cannot 
rule out the possibility that equation (Q) overestimates the value of R mq for an ADAF-magnetosphere 
interaction. The propeller effect would then be significantly less efficient than our previous results suggest 
(§5). 

If winds alone (p ~ 1) are to explain the reduced luminosity difference between BH and NS SXTs in 
quiescence, we need a highly inefficient propeller. For example, from the no-wind (p = 0) model to the 
p = 1 model described in §6.1, r mq increases by a factor ~ 5.8 which implies a~4x 10 4 times more efficient 
propeller according to equation ([ll]). Since the p = model has a propeller with / acc = 10 -3 , the propeller 
would have to be ~ 10 7 times less efficient in the p = 1 model than we would infer based on the theory 
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described in §5. In addition, in the presence of a wind, the value of / acc is likely only to be smaller than the 
expression in equation ([ll]). This is because the wind flows out preferentially from the polar regions of the 
ADAF (Narayan & Yi 1995a), so there is no gas to flow down to the polar caps of the NS. Even allowing 
for reasonable uncertainties in the scalings in equation (Q), it does not seem likely that the propeller model 
could be that wrong. We therefore feel that the p ~ 1 strong wind model is unlikely. 



7. Discussion 

Chen et al. (1998; their Fig. 3) argued that the segregation in luminosity swing between outburst and 
quiescence previously found by Narayan et al. (1997b) and Garcia et al. (1998) disappears in a large sample. 
We point out in this paper that the Chen et al. sample is larger only because they include numerous upper 
limits, many of which are based on observations with insufficient sensitivity to detect typical quiescent 
SXTs. Unless treated with care, for instance with a "detection and bounds" type method (e.g. Avni et al. 
1980; Schmitt 1985), upper limits can obscure a true correlation, and indeed this seems to be the case with 
the analysis of Chen et al. (1998). 

Our Fig. H shows a carefully selected and smaller sample, and we find that there is a difference of 
luminosity between quiescent BH and NS SXTs. Furthermore, we argue, based on theoretical mass transfer 
rates, that it is important to compare quiescent luminosities of NS SXTs and BH SXTs of similar orbital 
periods. We find that quiescent BH SXTs with short orbital periods are about 2 orders of magnitude fainter 
(in Eddington units) than quiescent NS SXTs with comparable orbital periods. The difference reduces 
to a little more than one order of magnitude when absolute luminosities are considered. The sample is, 
however, small. Future observations with Chandra X-ray Observatory (CXO) and XMM will hopefully 
increase the sample size sufficiently to provide a more convincing case. Hopefully, they will also allow 
separate comparisons of long and short period systems. This should help solidify the arguments advanced 
by Narayan et al. (1997b) and Garcia et al. (1998) for event horizons in black holes. 

We find that the luminosities of quiescent NS SXTs are lower than the values predicted by standard 
binary evolution models when combined with ADAF models. Following Zhang et al. (1998), we argue that 
the propeller effect offers a plausible explanation for the discrepancy. However, although our model for the 
propeller effect accounts for the luminosities and small emitting areas of quiescent NS SXTs, it suffers from 
several sources of uncertainty that may affect the results significantly. 

The predictions of the model strongly depend on the spin frequency f2* of the NS and therefore on the 
prescription chosen for the torque (Eq. |]l5| ). The magnitude of the torque is not well known, especially in 
the propeller regime (e.g. Davies, Fabian & Pringle 1979; Henrichs 1983). The results of our calculations 
are also sensitive to the precise values of the magnetospheric radii in quiescence and outburst and to the 
assumption that the fastness parameter w cr ;t = 1. 

If w cr it < 1, ri*, cq is reduced by a factor w cr it (Eq. [17]) and more mass is accreted because / acc is 
oc ui~f t . Note, however, that independently of equation [15| and the value of w cr ;t, spin periods of a few 
milliseconds and low magnetic field strengths (assuming that Eq. j?]] is correct) lead to predictions for 
luminosities and surface areas that are consistent with observations. 

Even if the spin periods of neutron stars in NS SXTs are known, the predictions of our propeller 
model are still somewhat uncertain because of the strong dependence of the fractions / acc and / SU rf on the 
magnetospheric radius in quiescence i? mq . If the standard expression (Eq. Q) overestimates the value of 
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-Rmq, then the propeller effect could be less efficient than indicated by our results in §5. In that case, winds 
from ADAFs could provide the additional mechanism required to explain the low quiescent luminosities of 
NS SXTs. 

We have also neglected the likely misalignment between the neutron star spin axis and magnetic axis 
and have assumed a dipolar structure for the magnetic field. Relaxing these assumptions could also modify 
our results, especially our estimates of the fractions / acc and / sur f- 

The system EXO 0748-676 appears unusually bright in our sample of quiescent NS SXTs. One possible 
explanation is that magnetic braking acts more efficiently in this system because of its rather short P rb, 
giving a larger mass transfer rate. If the mass transfer rate is higher, the propeller effect would also be less 
efficient, and would increase the mass falling on the NS even further. Note that EXO 0748-676 is known to 
be an unusual binary system, as shown by its unexplained variations of P or b (Hertz, Wood & Cominsky 
1997). Another possible explanation is that, for some reason (youth?), this particular system has not yet 
reached equilibrium and spins at a relatively slower rate than other NS SXTs. This would result in a less 
efficient propeller and a more luminous system in quiescence. 

This last point is a general prediction of our propeller model: the propeller effect is less efficient if the 
NS in a SXT spins slowly (P S pin ^ a few ms). It is possible that several such NS SXTs with slowly spinning 
NSs exist in our Galaxy. In these NS SXTs, most of the mass accreted via the ADAF would reach the NS in 
quiescence. The outbursts of these systems would therefore be of small amplitude, 2-3 orders of magnitude 
in X-rays rather than 5-6 orders of magnitude, because of their relatively high quiescent luminosities. 

In our model, the gas which is stopped by the propeller is neglected and its fate is left unspecified. 
This gas could affect both the dynamics and the emission properties of the accretion flow. For instance, it 
is unclear how the mass accreted via the ADAF will find its way out of the quasi-spherical accretion flow 
after being propelled outward. 

In §6, we considered ADAF models with winds assuming that the efficiency p of the wind is the same 
for BH and NS SXTs . In principle, p could be different for ADAFs around BHs and around NSs (because 
of the presence of the magnetosphere in NSs for instance) . Having a different value of p for BH and NS 
SXTs would not, however, affect our conclusion that winds alone probably cannot explain the observed 
quiescent luminosities of NS SXTs. 

Recently, Brown et al. (1998) argued that nuclear reactions in the crust of neutron stars, triggered 
during outbursts in NS SXTs, could efficiently heat up the NS cores in these systems. The energy deposited 
during an outburst would be reemitted during quiescence at a rate sufficient to explain the observed 
quiescent thermal emission. If so, the propeller effect must be even more efficient than our estimates in §5 
indicate, so that accretion in quiescence does not contribute much to the quiescent emission of NS SXTs. 

We note two criteria that bear on the relative importance of crustal heating and accretion to the 
quiescent luminosity of NS SXTs. First, in the scenario of Brown et al., the entire surface of the NS 
contributes to the quiescent emission. Small emitting areas have usually been inferred from the observations, 
by assuming that the quiescent emission is close to blackbody emission. However, Rutledge et al. (1998) 
have shown that much larger areas are inferred if detailed hydrogen atmosphere models are used instead 
of idealized blackbody models. There is a large uncertainty in these results because of the poorly known 
hydrogen column density to the systems. Nevertheless, the values of the emitting areas inferred by Rutledge 
et al for several systems are consistent with the entire NS surface being responsible for the emission in 
quiescence. This success favors the crustal heating model. 
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On the other hand, any observed rapid variability (e.g. on timescales of a few days) should be 
considered as strong evidence of the importance of accretion in quiescence, because the intrinsic thermal 
emission from the NS is not expected to vary on such small timescales (Brown et al. 1998). Campana et 
al. (1997) reported that the quiescent X-ray luminosity of the prototypical NS SXT Cen X-4 varied by a 
factor of ~ 3 over a period of 4 days. This establishes that the accretion luminosity can be more than 3 
times larger than any steady quiescent luminosity possibly present in this system. Though very little X-ray 
variability data are presently available, several studies of the optical counterpart of Cen X-4 in quiescence 
have shown a variability of order several tenths of a magnitude (e.g. Cowley et al. 1988; Chevalier et 
al. 1989; McClintock & Remillard 1990). Since it is natural to relate this optical variability to the X-ray 
variability observed by Campana et al. (1997), both sets of observations appear consistent with accretion 
being the dominant mechanism for the quiescent luminosity of NS SXTs. 

Throughout this study, we assumed that accretion is responsible for the quiescent emission of both 
NS and BH SXTs. The possibility that this emission could be due to two different processes in NS SXTs 
certainly complicates the comparisons between the two classes of systems. We note, however, that even 
if the dominant contribution to the quiescent emission of NS SXTs is crustal heating, the reason why 
quiescent NS SXTs, as a class, are brighter than quiescent BH SXTs remains the presence of a hard surface 
in NS SXTs and of an event horizon in BH SXTs. 



8. Conclusion 

In this paper, we have reconsidered the luminosity difference between black hole and neutron star soft 
X-ray transients in quiescence, which has been used to argue for the presence of event horizons in black 
holes. 

We show that the current observational data suggest that quiescent BH SXTs as a class are fainter 
than NS SXTs. This result agrees with the previous work of Narayan et al. (1997b) and Garcia et al. 
(1998), but it disagrees with the conclusions of Chen et al. (1998) who used an inappropriate sample of 
SXTs for their comparison. 

We point out that, for a reliable comparison of the two classes of SXTs, objects with similar orbital 
periods P or b should be compared. Otherwise, variations of the mass transfer rate Mt with P or b can mask 
the results. 

We find that the observed luminosities of quiescent BH SXTs are consistent with the predictions of 
binary evolution models for Mt if roughly a third of the mass supplied by the secondary is accreted by 
the black hole via an ADAF. This estimate is for oadaf = 0.3. The fraction goes down to 1/5 and 1/10 
for aADAF =0.1 and 0.025 respectively. The observed luminosities of quiescent NS SXTs suggest, on the 
other hand, that only a very small fraction of the mass transferred by the secondary reaches the neutron 
star surface. 

We explain the small fraction in NS SXTs by invoking an efficient propeller. We have constructed a 
model for the propeller effect that accounts for the observed luminosities of quiescent NS SXTs and their 
small emitting areas. In addition, the model appears to be consistent with the millisecond spin periods 
recently inferred from observations of accreting neutron stars in transient low mass binary systems (Aql 
X-l, SAX J1808.4-3658). 

Winds from ADAFs constitute an alternative explanation (BB98) for the very small fraction of mass 
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reaching the neutron star surface in quiescence. We argue that an ADAF model with strong winds and 
no propeller cannot explain the observed luminosities. However, an ADAF model with a wind of low or 
intermediate strength and a somewhat less efficient propeller (but still within the range of uncertainties 
of our propeller model) is consistent with the observed quiescent luminosities of BH and NS SXTs. In 
this case, the parameter S, which measures the fraction of the viscous energy that goes directly into the 
electrons, has to be large ~ 0.3 (Quataert & Narayan 1998). 

After submitting our paper for publication, we became aware of an independent paper by Asai et al. 
(1998) on the same subject. The conclusions of Asai et al. are very similar to those presented here: (1) 
observational data show evidence for a difference of X-ray luminosity between quiescent BH and NS SXTs, 
(2) only 0.1 — 1% of the mass transferred by the secondary reaches the NS surface in quiescent NS SXTs, 
and (3) this small accreted fraction could be due to the propeller effect in NS SXTs. 

The quiescent fluxes and upper limits reported herein are generally in very good agreement with those 
listed in Asai et al. (1998), with the exception of 4U1543-47. Here we assume a distance of 8 kpc, while 
Asai et al. assume 1 kpc. This accounts for most, but not all of the discrepancy between the Asai et al. 
upper limit and that in Orosz et al. 1998. We choose to use the latter limit, as it is more conservative. 
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Fig. 1. — (a) Quiescent luminosities L m ; n (Eddington units) in the 0.5-10 keV band of the NS SXTs (open 
circles) and BH SXTs (dots) listed in Table 1. The luminosities are plotted as a function of the orbital 
period P or b of each system, (b) Same as (a) except that L n ■ ^ a ci-,™™ in nni+o in38 
as (a) except that the ratio L m i D /L max is shown. 
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Fig. 2. — (a) Same as Fig la, except that only SXTs with well determined L m i n and P or b are shown, (b) 
Same as (a) except that L m i n are shown in units of 10 38 erg s _1 . (c) Same as (a) except that the ratio 

^min/imax IS shown. 
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Fig. 3. — Predictions of binary evolution models for the mass transfer rate Mr (in units of 10 18 gs" 1 ) in 
NS and BH SXTs, as a function of the orbital period P or b- GR refers to a phase of mass transfer driven by 
gravitational radiation and EXP to mass transfer driven by secondary expansion. 
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Fig. 4. — (a) The upper band shows the quiescent luminosities of BH SXTs in the 0.5-10 keV band predicted 
by ADAF models if all the mass transferred by the secondary is accreted via the ADAF. The lower band 
corresponds to ~ 1/3 of the mass transferred being accreted via the ADAF. This model fits the observed 
luminosities reasonably well. Both bands have been calculated for ccadaf = 0.3. The effect of varying «adaf 
is explained in the text, (b) The upper band shows the quiescent luminosities of NS SXTs predicted in the 
0.5-10 keV band if 1/3 of the mass transferred by the secondary reaches the NS surface. The lower band 
shows that the luminosities actually observed correspond to a very small fraction (<~ 10~ 3 ) of the transferred 
mass reaching the NS surface. 
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Fig. 5. — Fraction of mass transferred by the secondary that reaches the neutron star surface in the propeller 
regime, as a function of g = Madaf / M accura . The solid line shows this fraction for a neutron star with a 
surface magnetic field strength of 10 s G, while the dashed and dotted lines correspond to 10 9 G and 10 10 G, 
respectively. 
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Fig. 6. — (a) Fraction / sur f of the neutron star surface emitting radiation in the propeller regime, as a 
function of g — MADAF/A^accum- The solid line shows / su rf for a neutron star with a surface magnetic field 
strength of 10 8 G, while the dashed and dotted lines correspond to 10 9 G and 10 10 G, respectively, (b) 
Equilibrium spin period P sp in as a function of g = Madaf /Maccum for the same magnetic field strengths as 
in (a). 
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Fig. 7. — (a) The band shows the quiescent luminosities of BH SXTs in the 0.5-10 keV band predicted by 
ADAF models including a wind (p = 0.4, r tr = 10 4 , 5 = 0.3) if 1/3 of the mass transferred by the secondary 
flows into the ADAF. (b) The band shows the quiescent luminosities of NS SXTs predicted in the 0.5-10 
keV band if (1) 1/3 of the mass transferred by the secondary flows into the ADAF, (2) there is a modest 
wind (p = 0.4, r tr = 10 4 , S = 0.3) and (3) all the mass reaching the NS magnetosphere is accreted onto 
the surface. The inner radius of the ADAF (= the magnetospheric radius) is at r mq = 20 (see text). The 
wind reduces, but only by ~ one order of magnitude, the accretion rate at the magnetosphere (compare with 
Fig. 4b). We need a modest propeller to decrease the luminosity further to the observed levels. 
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Table 1: NEUTRON STAR AND BLACK HOLE SXTs 



System 


Porb (hi) 


D (kpc) 


log[i m in] (erg s *) 


mi(M ) 


(1) 


(2) 


(3) 


(4) 


(5) 


o EXO 0748-676 


3.82 a 


10 e 


34.1 


1.4* 


o 4U 2129+47 


5.2 fc 


6.3 e 


32.8 


1.4* 


o 1456-32 (Ccn X-4) 


15. l a 


1.2 e 


32.4 


1.4* 


o 1908+005 (Aql X-l) 


19 a 


2.5 e 


32.6 


1.4* 


o H1608-52 


98.4 C or 5 d 


3.6 e 


33.3 


1.4* 


• GRO J0422+32 (XN Per 92) 


5.1 


2.6 9 


< 31.6 


12 fc 


• A0620-00 (XN Mon 75) 


7.8 


l h 


31.0 


6.1 h 


• GS2000+25 (XN Vul 88) 


8.3 


2.7 h 


< 32.3 


8.5' 


• GS 1124-683 (XN Mus 91) 


10.4 


5 i 


< 32.4 


& 


• H1705-250 (XN Oph 77) 


12.5 


8.6' 1 


< 33.0 


4.9' 


• 4U 1543-47 


27.0 


8 e 


< 33.3™ 




• J1655-40 (XN Sco 94) 


62.9 


3.2^ 


32.4 


V 


• GS2023+338 (V404 Cyg) 


155.3 


3.5 h 


33.2 


12 h 



NOTE. - (1) o indicates a NS primary and • a BH primary. (2) Orbital periods from McClintock (1998), 
except where indicated. (3) Distances to the systems. (4) Luminosities in quiescence in the 0.5-10 keV band 
(corrected for the revised distances) from Narayan et al. (1997b) and Garcia et al. (1998), except where 
indicated. Quiescent luminosities for 4U1543-47 and GRO J1655-40 are based on 20.4 ks and 100 ks ASCA 
observations (respectively), and the quiescent luminosity for GRO J0422+32 is based on a 19 ks ROSAT 
observation. (5) Primary masses, (a) Van Paradijs (1995). (b) Simbad CDS Catalog, (c) Rittcr & Kolb 
(1998). (d) Chen et al. (1998). (e) Garcia et al. (1998). (g) Esin et al. (1998). (h) Narayan et al. (1997b). 
(i) Esin et al. (1997). (j) Hameury et al. (1997). (k) Bcckman et al. (1997). (1) Chen et al. (1997). (m) 
This is an arbitrary choice in the range 2.9 — 7.5 given by Orosz et al. (1998; see also Bailyn et al. 1997). 
(n) Orosz et al. (1998). (*) For simplicity, all NS masses are assumed to be 1.4 M Q . 



